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1 Abstract 

In this paper, we consider two basic questions about presenting a homoge- 
neous polynomial /: how many variables are needed for presenting /? How 
can one find a presentation of / involving as few variables as possible? We 
give a complete answer to both questions, determining the minimal number 
of variables needed, iVess(/), and describing these variables through their lin- 
ear span, EssVar(/). Our results give rise to effective algorithms which we 
implemented in the computer algebra system CoCoA | CoC04| . 

2 Introduction 

Polynomials, also seen as symmetric tensors, are ubiquitous in Applied Math- 
ematics. They appear in Mechanics ( |MQ02|), Sig nal and Image Processing 
f | CM96| ). Algebraic Complexity Theorv f |l3CS97| '). Coding and Information 
Theory r |Rom92j ). etc.. 

One of the main open issue is to manipulate polynomials in order to obtain 
presentations suiting the special needs of the application at hand. 

In Mechanics, it is often useful to separate variables. Given a polynomial 
f{xi, . . . , Xn), one splits the set of variables in two pieces, e.g. {xi, . . . , Xr} 
and {xr+i, . . . , Xn}, and a presentation of / is searched of the following type 

/(xi, . . . ,a;„) = g{xi, ...,Xr) + h{xr+i, ■ ■ .,Xn) 

for some polynomials g and h. 

Separating variables is a well established technique and the search for 
splitting methods in general is very active (see |MQ02| ). 

In Signal Processing, homogeneous polynomials (also known as quantics 
from ancient Invariant Theory) are of crucial importance. The main interest is 
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in the so called sum of powers presentations, where a homogeneous polynomial 
/ of degree d is presented as 

f = lt+... + lf 

where /i, . . . , are linear forms. 

Sum of powers presentations are treated in connection with quantics in 
| CM96| , while a more general approach relating them to Polynomial Interpo- 
lation and Waring Problem can be found in pIOl] . 

In this paper, we consider two basic questions about presenting a homo- 
geneous polynomial (from now on referred to as a form) in a "easier" way. 
Given a form /, how many variables are needed for presenting it? How can 
one find a presentation of / involving as few variables as possible? 

Even if these problems are so natural, we are not aware of a complete 
solution existing in the literature. In this paper, we give a complete an- 
swer to both questions. Our results give rise to effective algorithms which 
we implemented in the computer algebra system CoCoA (freely available at 
cocoa. dima.unige. it). 

More precisely, given a form / G S* = k[xi, . . . ,a;„], k any field, we call 
essential number of variables of f the smallest integer r for which there exists 
a set of linear forms {yi, . . . , yr} <Z S such that 

/ e k[yi, . . .,yr]; 

the linear forms yi, . . . ,yr are called essential variables of f. Then our main 
result is (see Definition^ Dcfinition|21and SectionOfor the notation involved): 

Proposition ^ Let f be a homogeneous element in S — k[xi, . . . , Xn] and 
T = fc[9i, . . . , dn] denote the ring of differential operators. Then 

Nessif) - rk{Cf), 

i. e. the number of essential variables of f is the rank of its first catalecticant 
matrix, and 

EssVar{f) = {Dof:DeTd-i), 
i.e. the essential variables of f span the space of its (d— 1)^^ partial derivatives. 

In Section 13 we briefiy recall some facts from Apolarity Theory which are 
the main tools of our analysis. 

In Section^ we use Apolarity and Catalecticant Matrices to obtain our 
main result. In Subsection 14. II we give some examples of the use of our algo- 
rithms. 

Remark 1. In this paper we work with forms, i.e. homogeneous polynomials. 
To apply our results to any polynomial /, it is enough to work with its ho- 
mogenization with respect to a new variable. Clearly, e.g., a presentation 
of f^ in essential variables readily produces a presentation of / in essential 
variables: it is enough to dehomogenize. 
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Remark 2. Throughout the paper k will denote a field of characteristic 0. 
Our results also hold in positive characteristic, but more advanced techniques 
are required for proving them, e.g. differentiations have to be replaced with 
contractions and divided powers have to be introduced (see |IK99n . 

The author wishes to thank B. Reznick and C. Ciliberto for their ideas 
on the problem. The CoCoA Team in Genoa, and especially Anna Bigatti, 
were of great help in the implementation of the algorithms. The comments 
and criticisms of the two anonymous referees were of help in improving the 
presentation of the results. 

3 Apolarity 

In this section we will briefly recall some basic facts from Apolarity Theory 
or, in modern terms. Inverse Systems Theory. Comprehensive references are 
|Ger96j . mainly Lecture 2, 6 and 8, and |IK99j . 
Consider the polynomial rings 

S = k[xi, ■■■,Xn] and T = fc[5i, . . . ,9„], 

where k denotes a field of characteristic 0, and give S a T- module structure 
via differentiation, i.e. we will think of T as the ring of differential operators 
acting on S. We denote this action with "o", e.g. dj o f = for f £ S. 

There is a natural perfect paring between homogeneous pieces of the same 
degree of S and T, namely 

Si X ^ k 

f D Dof^ 

is a perfect pairing for all «; in particular, St and Ti are dual to each other. 
Given subspaces, 

V CS.smdW C T, 

we denote by 

C and C Sr 

their orthogonal with respect to this pairing; notice, e.g., that dim^ y + 
dim/c dinifc Si = dim^ Ti. 

Given a form f E Sd, the ideal 

f^^{D€T:Dof = 0} 

is a homogeneous ideal of T and it is called the orthogonal ideal of /. 

Orthogonal ideals play a central role in the theory: they contain all the 
differential operators annihilating a given form and even more information, 
as it is shown by the following Lemma (for a proof see |Ger96j . Proposition 
8.10). 
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Lemma 1. Let f be a degree d form in S , then D ^ Ti,i < d, is such that 

Dof^O 

if and only if 

Do{D' o f) = 

for all D' £ Td-i- In other terms, for < i < d, {f^)i is orthogonal to the 
k-vector space spanned by the (d — i)^^ partial derivatives of f . 

Orthogonal ideals can be easily described introducing ad hoc matrices. 
In this paper it will be enough to describe the degree one part of a given 
orthogonal ideal, but similar descriptions exist in each degree. 

Definition 1. Let f d Sd and fix the standard monomial basis, e.g. with re- 
spect to lex order, {Mi, . . . , Mpf} of the k-vector space Sd-i. For i = 1, . . . ,n, 
consider the first partials 

dio f ^ CiiMi + . . . + CiNM^. 

The first catalecticant matrix of / is 

i = l,...,n,j = 1,...,N. 

For a general treatment of catalecticant matrices and their applications see 

iniMI, and HEM). 

Example L Let / — X1X2X3 6 k[xi, X2, X3] and consider the monomial basis 

{ xl,XiX2, Xi X3,xl, X2 X3,xl } 

of the space of degree two forms. Then 

/O 1 0\ 
= 1 . 
\0 1 0/ 

Catalecticant matrices determine the degree one part of orthogonal ideals 
readily: 

Lemma 2. Let f £ k[xi, . . . , Xn] be a form, then 

{aidi + . . . + a„a„) o / = 

if and only if the vector (ai, . . . , a„) is in the left kernel of Cf. In particular, 
dimfc(/^)i = n - rk{Cf). 
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Proof. The statement simply follows writing down the action of ai3i + . . . + 
a„9„ on / componentwise and considering the corresponding linear system of 
equations. 

Remark 3. Let / G . . . , a;„] be a linear form and consider its d-th power 
/ = l'^. Let L be a linear differential operator and notice that L o / = if 
and only if L o / = which is a linear equation in the coefficients of L. Hence 
rk(C;d) = 1 (actually, even the converse is true). In particular, this means that 
the form of Example ^ is not a pure power. 

4 How many variables? 

In this section we will use apolarity to answer our two basic questions: how 
many variables do we need to present a given form? How can we find a pre- 
sentation involving as few variables as possible? 

In what follows, we will work with the polynomial ring S = k[xi, . . . ^ x„], 
where k is any field such that char(fc) = (in positive characteristic similar 
results hold, but, in this paper, we decided to avoid the technical difficulties 
involved) . 

Lets introduce some definitions: 

Definition 2. Given a form f in S , the number of essential variables of /, 
Nessif), is the smallest integer r such that there exist linear forms yi, . . . ,yr (z 

5 for which f £ k[yi, . . . ,yr]- We call essential variables of / any set of 
generators of the k-vector space EssVar{f) = (j/i, . . . ,yr)- 

Roughly speaking, given a form / e 5, A'ess(/) tells us how many variables 
are necessary for presenting /, while EssVar(/) tells us how we can find such 
variables. In particular, it is clear that, if 

Nessif) = r and EssVar(/) = (yi, . . . , y^), 

then there exists 5 G . . . , j/^] C S such that f — g. 

Example 2. Consider the form / — /(xi, 0:2, X3) = [xi + X2){xi — x^)'^ in 
k[xi,X2, X3]. Clearly / is an element of the subring 2/2], where yi — X1+X2 
and y2 = xi — X3. Hence Nessif) ^ 2 and equality holds by Remark 13 
as rk(C/) ^ 1 and / is not a pure power. Also, notice that EssVar(/) — 
{xi+X2, X1—X3) and sets of possible essential variables are: {xi -I-X2, a;i — X3}, 
{x2 + X3, 2x1 + X2 - X3}, etc.. 

Using apolarity we can effectively determine A^ess and EssVar for a given 
form: 

Proposition 1. Let f be a homogeneous element in S — fc[a;i, . . . , x„] and 
T — k[di, . . . , dn] denote the ring of differential operators. Then 
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Nessif) = rk{Cf), 

i. e. the number of essential variables of f is the rank of its first catalecticant 
matrix, and 

EssVarif) = {Do f : D eTd-i), 
i.e. the essential variables of f span the space of its (d—V)^^ partial derivatives. 

Proof. If Nessif) then / e k[yi, . . . ,yr] for some linear forms yi,...,yr 
in S. Let 

(yi, ■ • -.Vr)^ = {Ll, ■ . ■,Ln-r) C Ti 

and notice that (/^)i 3 {Li, . . . , Ln-r)- Thus, by Lemma|21 we have rk(C/) < 
iVess(/). 

If rk(C/) = t, then = (Di, . . . , Dn-t)- Complete this to a basis of Ti 

Fi,..., Ft) 

and consider the dual basis of defined by the apolarity perfect pairing 

(zi, . . . ,z„_t,yi, . . .,yt). 

Hence, after a linear change of variables, we have / — f{zi, . . . , yi, . . . ,yt). 
But Dj annihilates all the elements of the chosen basis of 5*1 but Zj. As 
= (Z?!, . . . , Dn-t) we conclude that 

/ e k[yi, ...,yt] 

and rk(C/) > A^ess(/)- 

To conclude the proof, notice that the prefect pairing 5*1 x Ti ^ fc induces 
a well defined perfect pairing of fc- vector spaces 

where 

((/^)i)^ - (/:/e5i,Lo? = 0forallLe 

and, with the notations above, (^jr^ = O^i, • ■ ■ , Yt) and hence V — EssVar(/). 
The result follows applying Lemma = 1 case) which yields 

V={D'of:D'eTa-i). 

Example 3. Given the form 

f — x\ + X1X2 — 2xi X3 — 2xi X2 X3 + xi X3 + X2 X3 

we want to determine A'essC/) and EssVar(/). In order to apply Proposition 
n we compute the first catalecticant matrix of / 
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/ 3 2 -4 -2 l\ 
C/ = 1 -2 1. 
\-2 -2 2 2 0/ 

Hence -/Vess(/) — i'k(C/) = 2 and / can be presented as a form in two variables. 
To determine the essential variables of /, it is enough to compute the span of 
the second partial derivatives of /: 

EssVar(/) {x2 + X3,xi- X3). 

Summing these up, we see that there exists a degree 3 form g{yi, 1/2) G k[yi, 1/2] 
such that 

g{x2 +X3,Xi - X3) = f{xi,X2,X3), 

but how can we find gl 

To complete our analysis, we want to present a form / as a polynomial 
only involving essential variables: this can be done almost tautologically, but 
the notation are quite involved. We begin with an example. 

Example 4- Consider the form f ^ S ~ k[xi, X2, x^] in ExampleOl We already 
showed that there exists g G ^[2/1,2/2] C S such that f — g. To determine 
giui, 2/2), consider EssVar(/) — {x2 + X3, xi — X3) and complete its basis to a 
basis of 5*1: we choose {yi — X2 + X3,y2 — Xi — X3, zi = xi}. Hence we have 
a linear change of variables given by 

Xi = Zi, 

2^2 = yi + 2/2 - Zi, 
X3 = zi- 2/2. 

The basic requirement of the form (7(2/1, 2/2) is to satisfy the relation 

g{x2 +X3,Xi - X3) = f{xi,X2,X3). 

From this, changing variables, we get 

5(yi, 2/2) = f{zi, yi + y2- zi,zi- 2/2) = 2/12/2 + 2/2: 

which is the desired presentation in essential variables. As a byproduct, we 
readily see that 

f = {X2+ X3){xi - X3f + {Xi - X3)^ 

which is quite surprising considering the original presentation 

f = x\ + x\x2 — 2xf X3 — 2xi X2 X3 + X1X3 + X2 . 
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The procedure showed in the previous Example works in generaL Given a 
form / = /(a;i, . . . ,Xn) G S, we compute A^ess(/) = and we choose a basis 
for -/Vess(/) = (j/i; • ■ ■ j2/r); to avoid triviahty, assume r < n. Now, our goal 
is to determine g = g{yi, . . . , Ur) G ^[j/i, ■ • • , J/r] C iS such that f — g. To do 
this, complete the basis of Nessif) C 5i to a basis of 5*1 

Si = (j/i, . . . . . . ,Z„-r)- 

As Si — {xi, . . . ,Xn), the completed basis yields a linear change of variables 

(yi 7 ■ • ■ : Vr ; ^1 : ■ ■ ■ 7 ^n— r ) 7 
^niyi, ■ ■ ■ , Vr, Zi, . . . , Zn~r)- 

Notice that yi , . . . , are linear forms in S and hence there exist linear func- 
tions such that iji = yi{xi, . . . , x^), i = 1, . . . , r. Moreover, the following iden- 
tities hold by construction of (|) 

Vi = yi{xi{yi, ...,yr,Zi,..., Zn-r), ■ ■ ■,Xr{yi, . . . ,yr,Zi, . . . , Zn-r)) 

for i = 1, . . . , n. 

To determine g, it is enough to consider the desired relation 

/(xi, . . . , Xn) = 5(yi(a;i, • • ■ , Xr)-, yr{xi, Xr)). 

and to apply the linear change of variables (f). Thus we obtain g{yi, . . . , yr)- 

g{yi,...,yr) = 

= g{yi{xi{yi, ...,yr,Zi,..., Zn-r), ■ • • , Xr{yi, . . . , yr , Zi, . . . , Zn-r)), ■ ■ ■ 
. . . , yr{xi{yi,. ..,yr,Zi,..., Zn-r), Xr{yi, . . . , yr , Zi , . . . , Z^-r))) = 

— (yi ; • ■ ■ T yr, Zi, . . . , Zji—r): ■ ■ ■ i (yi i ■ ■ ■ ; yri ^1 ; • ■ ■ i Zn—r))- 

Notice that, as / and the functions Xi{yi, . . . , yr, zi, . . . , z„-r), i — 1, . . . , rt, 
are explicitly known, we have completely determined g as an element in 
k[yi, ■ ■ . ,yr]- 

Remark ^. As a straightforward application of the theory, we consider the 
detection of cylinders (i.e. algebraic surfaces ruled by a family of parallel lines 
moving along a fixed curve). Suppose you are given the polynomial equation 
of a surface J- : f{x, y, z) = in three space and you want to decide whether 
is a cylinder or not. It is well known that JF is a cylinder if and only if 
its defining equation is a function of two planes, i.e. there exist linear forms 
m{x, y, z) and l{x, y, z) such that /(x, y, z) — g{m, n) for some polynomial g. 
Hence, we readily have an effective procedure for cylinder detection: 

is a cylinder if and only if Aess(/'*) < 3, 

where denotes the homogenization of / (see Example [T)). Clearly, the 
method applies in any dimension for deciding whether a given hypersurface is 
a cylinder or not. 
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4.1 Using a computer 

The results of our analysis can be easily translated into algorithms and we 
wrote down procedures to be used with the Computer Algebra system Co- 
CoA. 

We begin with reporting a CoCoA session illustrating the use of our 
algorithms to work out the expository Examples |21 and 01 

Example 5. First we define the form we want to study 

F:=x"3 + x"2y - 2x~2z - 2xyz + xz~2 + yz"2; 

To compute the number of essential of variables, use the function NEssVar (F) : 

NEssVar(F) ; 
2 

To determine a choice of essential variables, use the function EssVar(F): 

EssVar(F) ; 

[y + z, X - z] 



Finally, NewPres (F) produces a presentation of the form involving the essen- 
tial variables y[l] = y + z, y[2] = x — z: 

NewPres (F) ; 
y[l]y[2]-2 + y[2]-3 



Usually, a given polynomial f{xi, . . . ,Xn) will essentially involve n vari- 
ables, i.e. -/Vess(/) = n. Hence our algorithms do not help in solving the 
polynomial equation f — 0. Nevertheless, our procedure should be used as a 
pre-processing tool. In fact, if the number of variables can be decreased, then 
the numerical solution of the equation can be performed much more efficiently. 
We illustrate this with the following "extreme" example. 

Example 6. We consider the degree three polynomial in four variables 

f{x,y,z,t) = fo{x,y,z,t) + fi{x,y,z,t) + f2{x,y,z,t) + fz{x,y,z,t), 

where 

/o = 3 

/i = -a; - y + 2^ + 3t 

/2 = 5a;2 4- lOxy + by^ - 20xz - 20yz + 20z'^ - SOxt - 30yt + 60zt + 45*^ 
f^^ + ix^y + 3xy^ + - Qx'^z - 12xyz - 6y'^z + I2xz'^ + 12yz^+ 

~8z^ - 9x'^t - ISxyt - 9yH + 36xzt + 36yzt - 36zH + 27xt^ + 27yf+ 

-hAzt^ - 27t^. 

In order to solve the equation f{x, y, z, t) — 0, we apply our algorithms to the 
degree 2 and 3 pieces of /: 
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EssVar(F2) ; 

[x + y - 2z - 3t] 



NewPres(F2) ; 
5y[l]-2 



and hence f2ix, y, z, t) — 5?/i^, where yi = x + y — 2z — 3t. Similarly 

EssVar(F3) ; 

[x + y - 2z - 3t] 



NewPres(F3) ; 
y[l]-3 



which yields fsix, y, z, t) = yi^. In conclusion, to solve the equation /(x, y, z, t) = 
0, it is enough to solve the equation in one variable 

yi^ + Syi^ - yi + 3 = 

and to apply some linear algebra to find all the solutions. 

We conclude with a Geometric example about cylinder detection. 

Example 7. Consider the degree five surface in three space J- : f{x,y,z) — 0, 
where 

/ = /O + /2 + /5 

and 

/o = -1, f2 = x^ -xy- 2y^ - iyz - z^, 

/g = a;^ + 2x^y - Ix^y^ - %x^y^ - Ixy'^ - 2y^ + ix'^z - \%x^y^z - 24:Xy^z+ 
-9y^z + Ix^z^ - Ylx^yz^ - SOxy^z'^ - IGy^z^ - Ix^z^ - \%xyz^+ 
-Uy'^z^ - ixz'^ - 6yz'^ - z^. 

In order to decide whether is a cylinder or not, we follow Remark ^ 
Introduce a new variable t and consider the homogenization of /, = 
t^fo + i^h + /s- Using CoCoA and denoting by FH the form f^{x, y, z, t), we 
get: 

NEssVar(FH) ; 
3 

EssVar(FH) ; 

[t, y + 2/3z, X + l/3z] 



NewPres(FH) ; 

-y[l]-5 - 2y[l]-3y[2]-2 - 2y[2]-5 - y [1] "Sy [2] y [3] - 
7y[2]-4y[3] + y[l]-3y[3]-2 - 8y [2] -3y [3] '2 - 2y [2] -2y [3] '3 + 
2y[2]y[3]-4 + y[3]-5 
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In conclusion, f^{x, y, z, t) = g{yi,y2, 2/3) where g is the output of the function 
NewPres(FH) and 

' yi=t 
y2 = y + lz . 
^y3 = x + \z 

Hence, we have the polynomial equality f{x, y, z) = g{l, y2, 2/3) and ^ is a 
cylinder ruled by lines parallel to the line 2/2 = 2/3 = 0- 
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